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I briefly review the formalism for treating the leading-twist two-gluon states ap- 
pearing in processes which involve r] and r)' mesons. The constraints on the size of 
the lowest order Gegenbauer coefficients of the two-gluon distribution amplitude 
are obtained from the fit to the r) and -q' transition form factor data. The results 
are applied to x — > wiv'v') decays and deeply virtual electroproduction of r] and 
7]' mesons. 



1. Introduction 

The description of tlie hard exclusive processes involving light mesons is 
based on the factorization of the short- and long-distance dynamics and 
the application of the pertm'bative QCD^. The former is represented 
by the process-dependent and perturbatively calculable elementary hard- 
scattering amplitude, in which meson is replaced by its (valent) Fock states, 
while the latter is described by the process-independent meson distribution 
amplitude (DA), which encodes the soft physics. The lowest Fock state of 
flavour-nonsinglet mesons consists of quark and antiquark, while for flavour- 
singlet meson the two-gluon state appears additionally. In this paper I give 
a status report on work of Ref. [2] , and discuss the proper treatment and 
the importance of these two-gluon states. 

On the basis of recent results'^'^ on rj and t]' mixing, we adopt the 
following representation of t] and rj' in a octet-singlet basis: 

\i] ) = cos 6*8 |?78) - sin 6*1 |r/i) , 



^On leave of absence from the Rudjer Boskovic Institute, Zagreb, Croatia. 
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\r]') = sin 6*8 |7?8) + cos 6*1 I771) , (1) 
where the pure octet ijjs) and singlet states are given by 

(t)i{x)\{uu-\- dd + ss)/\/^) + 4>g{x)\gg) , 

\m) = ^Mx)\{uu + dd-2ss)/V6), (2) 

and higher Fock states are neglected. In this representation, the mixing 

dependence is solely embedded in the and 0i angles, while in more general 
approach different distribution amplitudes (p^ and (p^ could be assumed for 
P = ri,r]'. The numerical values* /g = 1.26/^, /i = 1.17/^, Os = -21.2°, 
and 61 = —9.2° are used in this work. Alternatively, one could use the 
recently suggested quark-flavour basis'*, but the analysis of DA evolution is 
more straightforward in the above given octet-singlet basis. 




2. Two gluon distribution amplitude and the transition 
form factor for the flavour-singlet meson 

Employing for this section more transparent notation (pg = ^1, the DA 
evolution equation for rji takes the matrix form 

where ® denotes the usual convolution symbol. The kernel V is 2x2 matrix 
with a well defined expansion in as- The evolution of the flavour-singlet 
pseudoscalar meson distribution amplitude (DA) has been investigated in a 
number of papers^'^'''. Most of the results^'^ are in agreement up to differ- 
ences in conventions. On the other hand, the consistent set of conventions 
has to be used in calculation of both the hard-scattering and the distribu- 
tion amplitude, and these are not easy to extract from the literature. 

Following the recent analysis of the pion transition form factor®, we 
have performed a detailed next-to-leading order (NLO) analysis of the rji 
transition form factor taking into account both hard-scattering and per- 
turbatively calculable DA part, and this enabled us to fix and test the 
convention we are using, and to make a connection with other conventions. 

The hard-scattering amplitude one obtains by evaluating the 7* +7 — > qq 
and 7* +7 ^ 5<? amplitudes which we denote by Tqq{u, Q^) and Tgg(u, Q^), 
respectively. Owing to the fact that final state quarks and gluons are taken 
to be massless and onshell, Tqq and Tgg contain coUinear singularities, which 
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have to be factorized out in order to obtain the finite quantities and 

T{u,Q^) = (rH,g,-(x,g2,MF) TH,ggix,Q^,f4))(8Z-\x,u,f,l). (4) 

On the other hand, the unrcnormahzed quark and gluon distribution amph- 
tudes 4>q{u) and (f>g{u) arc defined in terms of (0|^(— z)7+75il4'(z)|r^i)and 
(0|G+"(— z) n G^+(z)|^i), respectively. The renormalization introduces the 
mixing of these quark and gluon composite operators and 

-^W = r L = -^(^' Mf) ® , _ • (5) 



We note that Z represents a 2x2 matrix. Perturbation theory cannot be 
used for a direct evaluation of $(w), but replacing \rji) by \qq) or \gg) enables 
us to obtain the perturbatively calculable DA part and to determine matrix 
Z. Finally, the rji transition form factor is given by 

Fj'fm iQ') = ^ T{u, Q2) ^ ^{u) = ^ Th{x, Q\fi%) ® ^{u, (6) 

Hence, the singularities, which appear in the calculation of the hard- 
scattering (4) and the DA part (5) should cancel, and we have used this 
requirement to check the consistency of our calculation. 

By differentiating (5) with respect to iJp one obtains the DA evo- 
lution equation (3) with evolution potential V expressed in terms of Z 
V = —Z~^ (g) (^Upd/diip Z). The solutions of the leading-order (LO) evo- 
lution equation are given by 



]{x,jj,%) = 6x{l - x) 



1 + B?,{u-l)C?,/'{2x-l) 



n=2 



,{X, ^,%) =X\1- Xf 'B'nil^F) CTli'^X - 1) , (7) 



n=2 



where 



2+ 

^niMp) - (MoJ ^ / 2 ^ ) + Pn (Mo) 1 ^ /, 2 ^ 
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The coefficients B^{fif,), i.e., Bn'"^\fif,), represent nonperturbative input 
at scale while 7± = 1/2 (7^' + 7^') ± ^ {7n - 7nf + ^ll^li^ and 



^59 1 ^19 
+ a > n — _ / n /Q\ 

^" ~ %+ _ ^59 ' ^" ~ 6 -7" - -v'* ^ ' 

in in " /n In 

are defined in terms of anomalous dimensions 7^^ = 7n^^, , 

The DA (/>g is normalized to 1, but, since l^dx 4'g{x,iJ,p) = 0, there is 
no such natural way to normalize 0g. It is important to emphasize that 
any change of the normalization of the gluon DA is accompanied by the 
corresponding change in the hard-scattering part. Namely, for (j)g — > cr^g, 
the projection of gg state on the rji state, which can be derived from the 
definition of the gluon distribution amplitude, gets modified by factor 1/cr, 
i.e., 

2 n- P x{l-x) (T 2 n- P x{l-x)' ^ ' 

and, thus, the hard-scattering amplitude changes while the physical quan- 
tity (F-y*-y^j , in this case) remains independent of the choice of convention. 
By inspecting Eqs. (7-9), it is easy to sec that the change of the normaliza- 
tion of (pg can be translated into the change of the off-diagonal anomalous 
dimensions 

7r-^7™ 7r-<T7r (12) 

and of the coefficient B~{^q) aB~{^Q). The former can be easily un- 
derstood in the "operator" language, i.e., by considering the impact of the 
change of the normalization of the gluon operator on the anomalous di- 
mensions. Hence, by employing (12) along with (10) and (11), we can 
consistently change our conventions in order to compare our results with 
other calculations from the literature. For historical reasons (comparison 
with the forward case), in what follows a = ^/nf/Cp is used. 



3. Applications 

Using the mixing scheme defined in Eq. (1), we have obtained the NLO 
leading-twist prediction for the r] and ij' transition form factors. For the 
treatment of ^s{x,lJ.p), we use the well-known LO result for the flavour- 
nonsinglet meson distribution amplitude^'^. We truncate the Gegenbauer 
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series at n = 2, and fit our results to the experimental data^. The fits are 
carried through with /i/? = and = Q, with as evaluated from the 

two-loop expression with n/ = 4 and = 305 MeV. For > 2 GeV^ 
and /io = 1 GeV, the results of the fits read 

^2(^*0) = -0-04 ±0.04 ^2^^^) = -0.08 ±0.04 = 9 ± 12 . (13) 

The existing experimental data and their quality allow us to obtain not 
more than a constraint on the value of i?| . As expected, we have observed 
a strong correlation between B\ and The quality of the fit is shown in 
Fig. 1. 
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Figure 1. -q (below) and 77' (above) transition form factors. The shaded area corresponds 
to the in Eq. (13) given range for B^{fi^) and BK/Jq). 



As a first application of the above given results, we have analyzed 
the Xco VV^ v'v' decays. Following previous work^*^, we take XcO as 
a non-relativistic cc bound state, and obtain the decay amplitudes for 
cc {qq){qq) and cc {gg){gg)- The ratio r,,'^'/r^^ is suitable for inves- 
tigating the sensitivity to the gluon contributions. Despite possibly large 
value of _B|, we have observed only modest dependence on the variations 
of Bi and Bg in the allowed range (13) (up to 20 % difference). 

Finally, I report on our result for the two gluon contribution to the 
deeply virtual electroproduction of 77, 77' mesons. In this case, the subprocess 
amplitudes q+Ji^ — *■ {qq)+q and g+7£ igg)+q have to be evaluated. Our 
result for the former amplitude H^^'^^^ is in agreement with the results from 
the literature^ In terms of subprocess Mandelstam variables (s,u,t — t), 
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the latter amplitude reads 



P(g) ^ 47ra, Cf 1 1 QV=^ T 0g(T) ^ /I 1_ 



In deeply virtual electroproduction of mesons (DVEM), the limit t — > 

has to be considered. For t = Q and s + u = — Q^, our result (14) gives 
-^"0+^0+ = 0. We conclude that the two-gluon contribution is suppressed. 



4. Conclusions 

We have performed a detailed analysis of the proper inclusion of the two 

gluon states in the (octet-singlet scheme based) description of the hard ex- 
clusive processes involving 77 and 77' mesons. Normalization and conventions 
have been fixed and discrepancies found in the literature have been resolved. 
From the fit of the r], rj' transition form factors to experimental data, we 
have obtained the range for the B2 , -Bf and coefficients of the (pi , (j>g 
and (j)s DAs, respectively. Expected strong correlation between B2 and Sf 
has been observed. The results have been applied to XcO VViv'v') decay, 
where only a moderate dependence on has been found. For DVEM, the 
1*Q ~^ {99)1 subprocess has been found to be suppressed for small momen- 
tum transfer t. We conclude that for considered processes the theoretical 
and experimental uncertainties do not allow further restriction of the pa- 
rameter range since only a modest dependence on the value of been 
observed. This contrasts the findings for rj' g*g* vertex^'^. 
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